Solar system tests of scalar field models with an exponential potential 
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We consider a scenario where the dynamics of a scalar field is ruled by an exponential potential, 
such as those arising from some quintessence-type models, and aim at obtaining phenomenological 
manifestations of this entity within our Solar System. To do so, we assume a perturbative regime, 
derive the perturbed Schwarzschild metric, and extract the relevant post-Newtonian parameters. 
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I. INTRODUCTION 

Scalar fields are a valuable tool in the development of 
proposals to address many of the fundamental questions 
in current physics, such as the ones associated to dark 
energy and dark matter. As is well known, the latter 
is required to account for the undetected mass disclosed 
by galactic rotation curves and cluster dynamics, while 
the former manifests itself through the accelerated ex- 
pansion of the Universe. In the context of quintessence 
models, most proposals require that the scalar field slow- 
rolls down an adequate, monotonical potential [1]; one 
way of implementing such behaviour relies on an expo- 
nential potential driving the dynamics of the scalar field. 
A scalar field with an exponential potential may also arise 
in the context of a dilatonic brane, due to its dynamics in 
the bulk [2]. Scalar fields are also central in alternative 
gravity models to account for the accelerated expansion 
of the Universe [3] and on the unification models of dark 
energy and dark matter [4, 5]. 

Even though scalar fields are primarily studied for their 
cosmological impact, they should also manifest them- 
selves on an astronomical scale and, more locally, within 
the Solar System. Although the derived effects are likely 
to be difficult to detect, constraints arising from local ob- 
servables may provide helpful insight into the inner work- 
ing of these scalar field models; one well known example 
is the Brans-Dicke gravity model [6] , where a scalar field 
plays the role of a variable gravitational coupling, and 
displays a non-canonical kinetic term. Current bounds 
(derived from the Cassini radiometric experiment) on the 
post-Newtonian parameter 7 (which equals unity in Gen- 
eral Relativity), I7- 1| < 2.3 x 10^^ [7], indicate that the 
Brans-Dicke coupling parameter obeys ui > 4.3 x 10^ (for 
a enlarged discussion of the latest bounds on General Rel- 
ativity and other fundamental symmetries, see e.g. Ref. 
[8]). 

Another example of the presence of an exponential po- 



tential lies in the so-called chameleon field model, where 
a scalar field is dynamically driven by a monotonical po- 
tential; in most applications of this model, an inverse 
power law potential, V{4>) = M'*+"0~", or fiducial po- 
tential V{(t)) = exp{Ml4)Y is assumed, with M a 
characteristic energy scale [9]; however, other examples 
consider that this potential assumes an exponential form, 
specially in the context of braneworld or quintessence de- 
rived chameleon fields [10]. The key feature of this model 
is a direct coupling with matter (through a conformally 
transformed "physical" metric), which enables a density- 
dependent effective potential possessing a minimum: this 
is akin to the Meissner effect in superconductors, with 
the non- vanishing minimum yielding a very small mass 
for the scalar field. Hence, deviations from expected be- 
haviour due to the chameleon field are usually posed in 
terms of this direct coupling with matter [9]. However, in 
this work we do not consider a matter distribution, but 
search for the vacuum solutions p = only (in this con- 
text, this is equivalent to neglecting the coupling between 
the chameleon and matter); for this reason, one falls back 
to the case of a dynamical scalar field driven solely by a 
monotonical exponential, and so we emphasize that the 
results do not apply to chameleon fields. 

In this work, we aim to establish the post-Newtonian 
effects arising due to the presence in the Solar System of 
a scalar field driven by an exponential potential, which 
assymptotically approaches a value given by its cosmolog- 
ical setting, defined as its profile when spacetime is char- 
acterised by a flat Friedmann-Robertson- Walker metric. 



II. THE MODEL 

One considers a model given by the Hilbert-Einstein 
action, containing a scalar field $ endowed with the usual 
kinetic plus potential terms, in addition to normal mat- 
ter. 
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S = I [kR^ + CJ\^f^ d'^x . (1) 
Jv 

where k, = (IGTrG)"-"^ and £$, Cm are, respectively, the 
Lagrangian densities for the scalar field $ and matter. 
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The former is given by 
1 



(2) 



where V{^) is the potential driving the dynamics of the 
scalar field. 

Assuming that temporal variations within the So- 
lar System occur at a timescale Hq^, one may take 
a quasi-static, spherical spacetime, with a Birkhoff-like 
Ansatz for the (non-isotropic) metric (with coordinates 
{t,r,9,(f>)), given by the lino clement 



-e^^df + a^t) (e^'^dr^+dn^) 



(3) 



where a is the scalar factor (with a{ta) = ao = 1), so that 
_ a^r^e^^'^ sinO] this is a simplifying assumption 
and neglects the more evolved issue of matching the inner 
Scharzschild-de Sitter metric with the outer Friedmann- 
Robertson- Walker (FRW), but will suffice for the purpose 
of this study, as shall be seen later. Clearly, the above 
metric asymptotically approaches the FRW metric if r — > 
and (7 ^ 0, as is expected. 

Also, notice that one cannot resort to Birkhoff's the- 
orem to justify the chosen Ansatz, since the presence of 
the scalar field indicates that one is not dealing with a 
vacuum field equation; indeed, the particular form for the 
metric is chosen as a natural candidate for the assumed 
quasi-static, spherically symmetric solution, neglecting 
non-diagonal terms in the metric. With this caveat in 
mind, the calculations below show that this solution is 
unique (up to second order, as shall be demonstrated), 
as it stems from a closed set of algebraic equations relat- 
ing the coefficients of the parameterised post- Newtonian 
(PPN) expansion of the functions r(r) and cr{r). 

Variation of the action with respect to the scalar field 
$ yields the Klein- Gordon equation = dV/d^; one 
begins by defining the "cosmological" profile of this scalar 
field, denoted by $c and defined by the Klein-Gordon 
equation on a flat FRW background. 



dV 



V! 



(4) 



where H = a/ais Hubble's expansion rate and dots de- 
note differentiation with respect to time. 

Given that the scalar field $ should approach asymp- 
totically its cosmological value one writes it in terms 
of a cosmological plus local contribution <[> = ^/'(r)+<I>c(t). 
With this decomposition, the D'Alembertian operator 
yields 



, (5) 



where the prime denotes diflFerentiation with respect to 
the radial coordinate r. Also, introducing an exponential 



potential V{^) = Vq exp{—l^), one obtains 

y($) = Vo e-'(*=+^) = Fee-'^ , 
with the definition Vc = Vqc"'*''; also, one gets 



dV 



Ve 



(6) 
(7) 



where = -ZVoe~'*= = -IVc- 

The Klein- Gordon equation then becomes 



-2a 



V'"+(^+r'-(7')V' 



(8) 

This is a generalization of a Yukawa type equation of 
motion in a curved spacetime, with a source term due 
to the absence of a defined minimum (that is, a non- 
vanishing V^'). At infinity, ijj ^ and the r.h.s. vanishes. 
Also, the usual Yukawa equation is recovered in the fiat 
spacetime limit. 

One now turns to the form of the relevant components 
of the Einstein tensor, namely 
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{2q - H)He-^'' + 
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(9) 



{2q-H)He 
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-2cr r 



t" + ij' - a') - + r' 



where q = —da/d^ is the deceleration parameter. 

The energy-momentum tensor of the scalar field $ is 



T 



SO that 



2 S{^U) _ 6U 



t; 



(10) 



(11) 
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V''^-F($) . 
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One assumes that normal matter behaves as a perfect 
fluid, with density p{r) and pressure p(r), related through 
an equation of state p = p{p). Its energy-momentum 
tensor is thus given by T/^ = diag{—p,p,p,p). One may 
now write the {t,t), (r, r) and (0,9) components of the 
Einstein equations: 
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- p-2<^ 2 
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T +(r-(7j -+T 

r 



= (47rG«>^ - Bil^^ + 47rG^i/>'2 + 87rGF($) + SttGp (12) 
= (^47rG$2 + (iJ _ 2g)ifj + 47rG-^V'^ - 87rGy($) + SttG^ (13) 

= e"^^ (47rG$^ + (i/ - 2g)if) - 47rG^^V''^ - 87rG\/($) + 87rGp . (14) 
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At infinity (with p(oo) = p(oo) = 0), these equations 
simpliiy to 

SttGYc = -47rG$2+3ij2 =4^Gr$2 + (if_2g)iJ . (15) 

For clarity, recall that the vanishing of pressure and den- 
sity refers only to the normal matter component; the the Einstein and Klein-Gordon equations, respectively: 



(positive) cosmological value of the scalar field's poten- 
tial Vc plays the role of a positive density pc and negative 
pressure Pc = —Pci as expected. 

Substituting Eq. (15), o(to) = ao = 1 and F($) = 

Vc exp{—lijj) into the above, one obtains the final form of 
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2^)+87rGp , 
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2")+87rGp , 


(b) 
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2^)+87rGp , 


(c) 








(d) 



(16) 



A fully analytic solution of the above system of differ- 
ential equations is not in general an easy task; for this 
reason, one resorts instead to an approximate solution, 
which should be self-consistent up to a desired order of 
accuracy. Since one aims at obtaining Solar System ob- 
servables which can be used to compare with predictions 
arising from the considered model, a natural candidate 
for a solution Ansatz is the so-called PPN expansion [7]. 
It should be noted that there is no guarantee that this 
is the only viable expansion of the full analytic solution; 
another possible approximation (not pursued here) could 
resort to a Yukawa-like metric, for example. 

The PPN expansion of the anisotropic metric 9^,^ and 
the scalar field ^ is written in terms of powers of U{r), 
the gravitational potential. The latter is defined as 



U{r, t) = G 



(17) 
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= 1 - 2J7 -I- 4a2t/ , 
= l + 2/3i?7 + 4/32f/2 
= l-F27i[/-F472f/2 



(18) 



up to second order in U (r) 
tion, this reads 



Outside the matter distribu- 
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l+/3r^ 
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-/32 - 



= 1-I-71 1-72 



with Rs = 2GM, so that 



(19) 



Outside the matter distribution, this collapses to the 

usual expression U{r) = GM/r, M being the central 
body's mass. The PPN expansion is given by 



1 

2 r 2 



2a2 



(20) 
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Rs 3(1-2^2 ) /^i?. 
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-71^ + (7?- 27.) 



271^+3(272-7?) f- 



In order to obtain a consistent solution outside the cen- 
tral body = /> = 0), one could proceed and substitute 
the above expansions into the set of Eqs. (16), equating 
the coefficients of same order terms in i?s/r, and solving 
the resulting set of algebraic equations for the parame- 
ters aj,/3i,7i. However, inspection of Eqs. (16) shows 
that the leading terms in the l.h.s. are of order or 
r~^; on the r./i.s., the ^/^'^ leading term is of order r~^. 
The second term on the r.h.s. may be expanded as 



-lib -It 
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(21) 
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Therefore, one may first ensure that no terms of order 
above are present, by solving the two algebraic equa- 
tions that follow from the vanishing of the coefficients of 
the above expansion. This yields 



71 



-1 



72 = 0:2 



(22) 



which, from the expansion Eqs. (20), amoimts to i/' = 
2t//; also, it follows from the above reasoning that this 
is the only approximate PPN solution for Eqs. (16), up 
to second order in Rsjr. 

Since the gravitational potential (and first derivative) 
is continuous across the central body's surface, so are the 
metric components 500 and g^r (for any set of PPN pa- 
rameters cii , l3i). The continuity of the scalar field is more 
troublesome; indeed, the above solution ip = 2t/1 is valid 
outside the matter distribution, but may break inside it, 
where terms dependent on the matter variables p and p 
play a relevant role. Also, accounting for a possible cou- 
pling between the scalar field and ordinary matter would 
amount to an even more involved interior solution for the 
scalar field. Hence, the desired continuity across the cen- 
tral body's surface is not regarded as a consistency test 
of the validity of the approach taken in this work, but of 
the eventual coupling to be included in the full theory. 
Regarding this point, it should be noted that, since both 
the density and pressure vanish at the surface, but the 
derivatives do not, a coupling function that depends only 
on p and p should be preferred, given the form of Eqs. 
(16). 



To evaluate the accuracy of this solution, one resorts 
to the Klein-Gordon Eq. (16-rf), which now reads 



yielding the constraint 

-l-H4a2 + /3i =0 , 



(23) 



(24) 



required for the solution ^ = 2r// to be valid up to sec- 
ond order in Rs/r. 

With this solution, the Einstein equations outside the 
central body become 

1 - e-2<^ 2 



-^'e-'^=^e--r- , (a) (25) 

-^-^ + = Ae-2V2 , (6) 

r" + {r'-a')(^l+r''^=-Ar" , (c) 

where the dimensionless constant A = WnG/l^ = 
167r(M/Mp;)2 has been introduced, with M = being 
the characteristic energy scale of the scalar field. Adding 
Eqs. (25- a) and (25-6), one gets 



-{T' + a') = AT'^ . 
r 

Adding this to Eq. (25- c) yields 

t"+(^+t'-c7'\t' = , 



(26) 



(27) 



which is just the Klein-Gordon equation (16-d) with the 
solution ip = 2t/1. Hence, one only retains Eqs (25- a, b), 
obtaining 



1 (=2(7 _ 1 
a' = tArr'^-- 

2 2r 

1 - 1 

t' = -ArT'^ + - 

2 2r 



(28) 



Introducing the metric coefficients Eqs. (20) and tak- 
ing terms up to second order in Rs/r yields the set of 

equations 

/3i R, (Pi \Rl 13, R. (A Rl 
- — ^ + I — - ^2 ) ^ - - + I ^ - — ) — > 



2 r2 V 2 
R. fl 



i?2 (A (32\ H] 



2 / r3 



Equating first and second order terms and solving the 
resulting set of algebraic equations yields 

/3i = l , /32 = l-j , a2=0 . (30) 
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These results satisfy the constraint Eq. (24). In the 
absence of the scalar field <C 1 ^ ~ 0), one re- 
covers the usual expansion of the Scharzschild metric in 



anisotropic coordinates, goo s 



-(1 - Rs/r) 



and grr s = e^"" = [l- Rs/r)~'^. Notice that the ef- 
fect of the scalar field is very dim, of second order in the 
radial component of the metric. 



III. PPN METRIC 

To read the PPN parameters, one must transform to 
an isotropic metric of the form ds^ = hoodt^ + 

hpp [dp^ + p^dfl^). Equating the line element to the 
anisotropic form ds"^ = goodt^ + grrdr"^ +r^dri^, one gets 
goo{r) = hoo{p), grrdr^ = hppdp^ and = hppp"^. Divid- 
ing these and taking the square root gives 



1/2 

Irr J "■r 

ar = — 

r p 



dp 



log p ■■ 



It 



dr + const. . (31) 



Since one expects a small perturbation to the stan- 
dard isotropic coordinates po (derived from the usual 
Scharzschild metric), it is advantageous to write (drop- 
ping the integration constant) 



log p 



dr 



\/r'^ — RsT 

Expanding the exponential to second order, one gets 



(32) 
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yielding 



r 



log p ■ 



1 



dr 



(33) 



(34) 



'•^ — RsT 

The standard radial isotropic coordinate, defined by 

(35) 



dpo 

PO 



dr 



y/r"^ — Rgr 



is given by 



Po 



Rs 
r 

so that 
log p = 




(36) 



y/r'^ — Rgr 



dr = 



(37) 



Po 8 pI 



log Po + 



A 



dpo 
2 



(1 + 4-)^ 3(1 + ^), 



log Po 



A fR 



16 \po 



enabling to rewrite 



p = Po exp 



A_ /R, 
16 \po 



(38) 



and hence, 



Po 



16 V P 



(39) 



Finally, one gets 



R 1-A /R, 
2p 16 \ p 



(40) 



One can now write the isotropic metric components: 

R« 



-hoo{p) 



-9oo{r{p)) = 1 



r{p) 



(41) 



and 



hppip) 



l-2C/ + 2i72-^i^[/3 



Rs 1-A f Rs 



2p 16 



l + 2U+^—^U^ 



(42) 



where U{p) — Rs/2p. Comparison with the standard 
PPN isotropic metric yields /3 = 7 = 1, indistinguishable 
from General Relativity. This is because the scalar field 
manifests itself only in second order in the Scharzschild 
metric radial component. However, non- vanishing effects 
are present beyond post-Newtonian order. 



IV. POTENTIAL OBSERVATIONAL EFFECTS 

One may discuss the magnitude of A that would still 
lead to a possibility of detection, by considering the cur- 
rent sensitivity of experiments; forfeiting a thorough dis- 
cussion of available techniques and instruments, we con- 
sider as figure of merit a e = 10~^^ m/s^ sensitivity for 
acceleration measurements (see, e.g. Refs. [11] and [12]). 
For instance, a relative frequency sensitivity of 10""'^^ is 
achievable; however, since = \/ goo{A) / goo{B) , 
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50.0 




r/lAU 



FIG. 1: Exponential potential mass scale M against mini- 
mum detection distance from the Sun vm, for an acceleration 



sensitivity e = 10 m/s 



these measurements are only sensible to changes in goo, 
which is unaltered up to second-order in i?,,/r. 

The radial acceleration (using the isotropic metric) is 
given, in the Newtonian limit, by 



00 



00 5 



(43) 



so that the contribution due to the presence of the scalar 
field $ is 



— — 



5 {GMf 



4 {rcy 



A = 1.5 X IQ-^A m/s^ 



(44) 



where one assumes an acceleration measurement very 
close to the Sun's surface, r ~ Rq « 7 x 10^ to. Equat- 
ing this with the mentioned sensitivity e, one concludes 
that the dimensionless parameter A may be detected 
within the range A > 6.5 x 10"''; using the defini- 
tion A = 16nG/P = 16TT{M/Mpi f, this implies that 
M > 3.6 X 10-^ Mpi fa 4 X IQis GeV/c^. Exponentially 
driven scalar fields arising in quintessence models usually 
display M Mpi , well within the obtained range. 



Inversely, one may obtain the relation between the 
mass scale of the exponential potential M and the maxi- 
mum distance vm at which acceleration measurements 
may detect the effect of this scalar field. From Eq. 
(44), one obtains tm = 7.8 x IQ-'^y/M/Mpi AU, or 
M/Mpi = 166(rM/l ^^7)^ as plotted in Fig. 1; for 
M Mpi, one gets about one fifth of Mercury's semi- 
major axis. 

Notice that, although this effect may be detectable, 
many issues common to space science experiments should 
be accounted for, namely the need to clearly discrimi- 
nate between several spurious contributions to the over- 
all acceleration of a test body - such as solar pressure, 
environmental effects (magnetic fields, drag from inter- 
planetary dust distributions, etc.), competing relativistic 
effects, and many others (see e.g. Ref. [12-14] for a 
discussion). 

V. CONCLUSIONS 

In this work we have addressed the issue of deriving 
the local effects of an exponentially driven scalar field. 
The latter could arise from a variety of sources, including 
quintessence models [1] , dilatonic branes [2] , Brans-Dicke 
models [6] or the chameleon field [9], amongst others. We 
have obtained results which are valid up to second order 
in the gravitational potential U; these indicate that no 
manifestation of such a scalar field should be visible at 
post-Newtonian order, but arise at post-post-Newtonian 
order. If the expected range of the exponential poten- 
tial is of the order of Mpi , we conclude that it may be 
possible (although experimentally challenging) to detect 
these feeble effects in the vicinity of the Sun. 
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